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A general expression for the steady-state matrix cracking stress in reinforced brittle matrix composites is derived using a 
/-integral analysis. The result is expressed in terms of the stress-displacement relation that characterizes the stretching of 
crack bridging ligaments. The influence of residual stress is assessed and a condition for spontaneous matrix cracking due to 
residual tensile stress in the matrix is- evaluated. Results of the analysis are compared with independent solutions for 
composites containing unbonded reinforcing fibers. 



1. Introduction 

If the strengths of crack bridging ligaments in a 
composite material exceed a critical value, a crack 
can be made to extend indefinitely in the matrix 
without the ligaments rupturing in its wake. Fur- 
ther loading causes formation of periodic matrix 
cracks, with separation dictated by a characteristic 
transfer length associated with the bridging liga- 
ments (Aveston and Kelly, 1973; Aveston et al M 
1971). This failure mechanism has been observed 
in several important composite systems, including 
reinforced cements, glasses, and glass ceramics 
(Aveston et al M 1971; Marshall and Evans, 1985; 
Sambell et ah, 1972; Phillips, 1972; Brennan and 
Prewo, 1982; Prewo and Brennan. 1982; Ali and 
Grimer, 1969; Majumdar, 1970; DeVekey and 
Majumdar 1968; Allen, 1971). An analogous 
mechanism has also been observed in connection 
with the cracking of thin glass films on metal 
substrates (Raj, 1987). Th e stress required to ex- 
tend the first matrix crack (i.e.. the onset or 
damage) is a decreasing function of the crack siz e 
jfjh? crack is sm^lL but app roaches a_ constant 
"steady-s tate" value JbHarge cracks (Marshall et 
al,T$85). In this paper we present a simple deriva- 
tion of a general solution for the steady-state 
matrix cracking stress. The solution, which is ob- 
tained by use of the /-integral (Rice, 1968), is 



equivalent to one obtained recently by Rose (1987), 
who evaluated changes in configurational energy 
involved in transporting a strip of material from a 
location far ahead of the crack tip to a position far 
behind the tip. 

The steady-state matrix cracking stress has been 
calculated previously for some specific, relatively 
straightforward ligament bridging mechanisms. 
Aveston, Cooper, and Kelly (1971) first analyzed 
bridging that involves frictional sliding between 
fibers and matrix, in the limit of small frictional 
stress. More recently. Budiansky, Hutchinson, and 
Evans (1986) obtained ji general solution fQ LiQ£i 
uonal sliding , valid over the full range of frictional 
stresses, from the small stress limit of Aveston, 
Cooper, and Kelly to the large stress limit at 
which sliding does not occur and the bridging 
ligaments act as linear springs. These analyses 
involved somewhat detailed calculation of energy 
changes, similar to the approach adopted by Rose 
(1987). An alternative, stress intensity approach 
has been used to evaluate both steady-state and 
nonsteady-state cracking for frictional and linear- 
spring ligaments (Marshall et al., 1985; Marshall 
and Cox, 1987). However, evaluation of the 
steady-state matrix cracking stress by this ap- 
proach entails numerical solution of an integral 
equation. An appealing feature of the solution 
discussed in this paper (in addition to its simple 
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analytic form and generality) is that it is expressed 
directly in terms of the stress-displacement rela- 
tion that describes the stretching of the bridging 
ligaments. This form of the solution allows changes 
in matrix cracking stress due to changes in the 
bridging mechanism or to the presence of residual 
^stresses to be readily deduced. 

2. General solution for steady-state matrix crack- 
ing 

The opening of a crack that is bridged bv 
reinforcements involves stretching of ligaments 
between the crack surfaces. This stretching may be 
characterized by a relation between the stress a, 
in the ligament and an average local crack 
displacement, u, as depicted in Fig. 1. The form of 
this relation depends on the details of the bridging 
mechanism and reflects properties such as liga- 
ment deformauon, reinforcement/matrix debond- 
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(c) 



F.g. 1. (a) Crack bridged by reinforcing ligaments; (b) Sche- 
mat.c of stress-displacement relation for stretching of bridging 
ligaments; (c) Crack with bridging ligaments replaced by 
surface tractions. y 



>ng and factional sliding, as well as elastic stretch- 
ing of the reinforcement. The peak value a t = S 
represents the strength of the ligaments. The ' 
decreasing portion of the a,{u) curve depends on 
the nature and location of reinforcement failure 
this is important for determining the toughness of 
composites in which reinforcements rupture in the 
wake of the crack, but i, does not influence 
steady-state matrix cracking. 

The influence of bridging ligaments on the stress 
anywhere in the body can be evaluated bv replac- 
ing the ligaments with crack surface tractions equal 
•n magnitude to the stress. o r . in the stretched 
hgament (F.g. 1(c)). To proceed further, the com- 
posite is approximated as a continuum with the 
average elastic properties of the composite, and 
with continuous pressure 

P(u)=fo t (u) (1) 

acting over the crack surfaces (/ is the volume 
traction of reinforcements). This continuum ap- 
proximation requires the crack to be large com- 
pared with the dimensions of the microstructural 
features that give rise to the bridging ligaments. 
We consider a crack with intact ligaments over its 
entire surface and subject to remotely applied 
uniform, normal, tensile stress. *„ that increases 
monotonically from zero to a value smaller than 
/i. The crack opening displacement and the crack 
surface pressure increase monotonically with dis- 
anoe behind the crack tip and, for sufficiently 
long cracks, approach asymptotic limits equal to 
u and a (p cannot exceed a,) at the mouth of 
he crack (Marshall et al„ 1985) (Fig. 2(a)). This is 
the steady-state matrix crack configuration. The 
stresses at the crack tip increase as the applied 
stress increases, but are independent of the total 
crack length. 

For the purpose of calculating the stresses and 
strains near the crack tip, the remotely applied 
uniform stress, •„ can be replaced by a uniform 
opening pressure acting over the crack surfaces. 
Then the resultant crack surface pressure becomes 
o = a i -p( u ) (2) 

as depicted in Fig. 2(b). This pressure is maximum 
at the crack tip and. for the steady-state matrix 
crack, approaches zero far from the crack tip 
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Fig. 2. (a) Steady-state matrix crack loaded with uniform 
applied stress; (b) Steady-state matrix crack with uniform 
applied stress replaced by uniform opening pressure acting on 
crack surfaces. 



The crack configuration of Fig. 2(b) is conveni- 
ent for evaluating the /-integral (Rice, 1968). For 
the closed path shown in Fig. 2(b), there are three 
contributions that must sum to zero: 

4+^ + ^,ip = 0. ( 3 ) 

The term / x from the path armmoLtJie^t ress-free 
outer bou ndaries of the body is zero. The contri- 
bution from the path over the"crack surfaces is 
(e.g. Lawn and Wilshaw, 1975) 



J B = ~2f\(u)du. 



(4) 



and the contribution from the path around the 
crack tip can be expressed 

'up = * 2 (W 2 )/£, ^, r (5) 

where K is the crack tip stress intensity factor and 
E and v are the Young's modulus and Poisson's 



ratio for the composite. Combination of eqns. 
(2)-(5) yields the result 

KHl-p 2 )/2E^o^ r p (u)du (6) 

Jq 

= C ud P- (7) 

Equation (6) is a general expression relating the 
crack tip stresses (or 7 up ) to the applied stress. It 
is noted that this expression is equivalent to the 
result derived by Rose (1987) and that the right- 
hand side is the complementary energy density. 

The critical applied stress, o c , for steady-state 
matrix cracking can be evaluated for any given 
bridging law by setting / lip =/ c (or equivalently 
K=Kc) as the criterion for crack growth in the 
matrix. However, it is necessary to express this 
fracture criterion in terms of the fracture tough- 
ness of the unreinforced matrix. Two approaches 
have been used for this. In one, the fracture en- 
ergy is taken to be reduced by the factor (1 - 0 to 
reflect the true area of fracture surface created in 
the matrix (Budiansky, 1986; McCartney, 1987: 
Kelly and McCartney, 1987), i.e. 



■Wc m (W) 



(8) 



where J c m is the fracture energy of the unrein- 
forced matrix. An alternative (upper bound) crite- 
rion was obtained by taking the stress intensity 
factor in the matrix to be equal to the critical 
stress intensity factor, K c m , in the unreinforced 
matrix (Marshall et al., 1985), i.e. 



K C = K?E/E U 



(9) 



where £ m is the Young's modulus of the matrix, 
and the factor E/E m accounts for the relative 
magnitudes of the average composite stresses and 
the actual stress in the matrix. With eqn. (5) this 
criterion can be expressed 



J,-E{Kn\\-^)/El 



(10) 



which differs from eqn. (8) by the factor £ m (l - 
/)/£ (since 7 c m = (K?) 2 (\ - v 2 )/E m ). The dif- 
ferences between these two criteria are discussed 
more fully elsewhere (Marshall et a!., 1988). For 
consistency with other publications, we shall adopt 
the criterion of eqn. (8) in the following sections. 
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3. Comparison with alternative solutions: Un- 
bonded fibers 

To evaluate the energy changes accompanying 
steady-state matrix cracking in unbonded fiber 
reinforced composites. Budiansky. Hutchinson and 
Evans (1986) calculated the stresses and strains in 
the bridging fibers and the surrounding matrix 
using a modified shear lag analysis. By integrating 
these strains to obtain the crack opening displace- 
ment, the stress-displacement relation for the 
bridging fibers was recently derived in the follow- 
ing form (Marshall et al.. 1988): 

■ (*' P<P*.u<u* 
\fi(p 2 + a*/4fi*) p>p* <u>u * (") 

where 
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P* = a/2B 

u* = cr/2B 

REJl-f) 
pfE t E 

*^(1 -ft 



P 2 = 



4E,E 2 / 2 t 



-4E(1-/) 



£ f(l + " m )[2 log/+(l -/)(3_/)j 



(12) 
03) 

(14) 

(15) 
(16) 



with R the fiber radius, E, the Young's modulus 
of the fibers, » m the Poisson's ratio of the matrix, 
and t the sliding frictional stress at the fiber 
matrix interface. The condition (u=u m , p=p*) 
defines the onset of sliding between the fibers and 
matrix. For displacements smaller than «*, sliding 
does not occur and the fiber bridges act as linear 
elastic springs. For u»u*, sliding between the 
fibers and matrix occurs over an area that extends 
a large distance (compared with R) from the crack 
surface and the solution corresponds to the large 
slip limit analyzed in the steady state by Aveston 
Cooper and Kelly (1971). 

3.1 No-slip limit: Linear springs 

If the matrix cracking stress is smaller than p* 
slip does not occur at any of the crack-bridging 



fibers. The stress-displacement relation is (eqn. 

" = ap (17) 
and substitution in eqn. (6) gives the solution for 
the steady-state matrix cracking stress 

(W 2 . (18) 



With J and a defined by eqn. (8) and (14). it is 
straightforward to show that eqn. (18) is identical 
to the energy balance solution of Budiansky 
Hutchinson and Evans (1986). 

3.2 Large slip solution 

U the matrix cracking stress is much larger than 
P . the bridging stresses are dominated by the 
large displacement limit 



u = Bp 2 . 



(19) 



In this case the steady-state matrix cracking stress 
obtained by substitution in eqn. (6) is 



»i-(3V20)' 



(20) 



With J c and B given by eqns. (8) and (15), this 
result is identical to the large slip solution of 
Aveston, Cooper and Kelly (1971). 

3.3 General solution 

The general solution for steady-state matrix 
cracking stresses, o c , larger than p* is obtained 
by substituting the complete expression of eqn. 
(11) into eqn. (6) yielding 

• /c = (¥)[°c 3 + 3p*a c -(^) 3 J. (21) 

This result can be expressed in terms of the limit- 
>ng so utions o 0 and a, for no slip and large slip 
by making use of the relations 

20o> o 0 3\a 0 l (22) 

which follow from eqns. (12), (18) and (20). The 
resulting expression, 
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(20) 



(21) 

t limit- 
ge slip 



. (23) 



is identical to the energy balance solution of 
Budiansky, Hutchinson and Evans (1986). 



4. The influence of residual stress 

Many composites contain residual stresses due 
to differences in nonelastic strains that occur in 
the matrix and reinforcements during fabrication 
(e.g., thermal expansion mismatch, plasticity, 
phase transformation). The derivation of eqn. (6) 
is based on a continuum model that requires all 
dimensions to be large compared with the micro- 
structural features that give rise to the residual 
stresses and the bridging ligaments. Over such 
dimensions the average residual stress is zero. 
Therefore, the derivation of eqn. (6) is not altered 
in any way by the presence of residual stresses 
(provided, of course, that further phase transfor- 
mation or plasticity does not accompany matrix 
cracking). However, residual stresses do influence 
a c through their effect on the stress-separation 
function p(u). 

The influence of residual stress on the relation 
p(u) has been evaluated for composites with un- 
iaxially aligned fibers (Marshall and Evans, 1988). 
In this case the longitudinal residual stresses in the 
matrix, , and fibers, a^, are related by 



(24) 



One effect of the residual stress is to shift the 
intercept of the function p(u) along the stress axis 

by 



a°=-o?£/£ m 



(25) 



as depicted in Fig. 3. In general, residual stress 
can also alter the slope of the p(u) relation. 
Therefore, this relation can be written in the form 



p(u) = ai + />,(«. Or), 



(26) 



where the subscript r denotes a modified function 
in the presence of residual stress, and eqn. (6) 
becomes 



K 2 (\-u 2 )/E = 2(o a -<,° K )u> 

-if > r («, a R )d« 



(27) 



og 1 COMPRESSIVE 



ZERO RESIDUAL 
STRESS 



P<u) 




Fig. 3. Influence of residual stress on the stress-displacement 
relation for ligament stretching. 



For certain bridging mechanisms, the effect of 
residual stress is to translate the function p(u) 
uniformly along the stress axis. Two examples are 
fibers that do not experience debonding or sliding 
and fibers held by friction due to surface rough- 
ness at the fiber matrix interface (i.e., the sliding 
resistance, t, independent of residual stress). In 
this case 



(28) 



where p 0 (u) is the stress-displacement function in 
the absence of residual stress, and both u a and the 
integral in eqn. (27) are independent of a£. There- 
fore, the effect of the residual stress is simply to 
increase the matrix cracking stress by o£ (a£ is 
positive for compressive residual stress in the ma- 
trix and negative for tensile residual stress). 

Equation (27) with K=K C provides a general 
solution for the critical applied stress for matrix 
cracking in the presence of residual stress. Alter- 
natively, we can solve for the critical residual stress 
that will cause spontaneous matrix cracking by 
setting o a - 0 at K - K c . (In this case o£ will be 
negative, i.e., tensile residual matrix stress.) 



5. Discussion 

An insightful representation of the general rela- 
tion for steady-state matrix cracking (eqn. (6)) is 
shown in Fig. 4(a). The right-hand side of eqn. (6) 
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STEADY-STATE 
MATRIX CRACKING 
STRESS 



AREA * A J c /2 

= STEADY-STATE 
TOUGHNESS 
INCREASE 




F* 4 Representation of (a) steady-state matrix cracking stress 
and (b) steady-state toughness increase, in terms of areas* 
related to the stress-displacement curve for ligament stretching. 



is given by the shaded area between the stress-dis- 
placement relation for the bridging ligaments and 
the constant stress line representing the applied 
stress. The critical condition for matrix cracking is 
determined by the applied stress for which this 
area is equal to JJ2 (i.e. K?(l - „ *)/2£). Thus, 
for a given matrix and volume fraction of rein- 
forcement this area is constant. It also follows that 
any change in the matrix toughness or composite 
stiffness that increases the ratio KJE must cause 
o e to increase. Moreover, the effect of changing 
the nature of the bridging ligaments on the matrix 
cracking stress can be readily deduced; generally 
any change that stiffens the loading portion of the 
curve p(u) must increase o c , whereas changes to 
the maximum value of p(u) or to the region of the 
curve beyond the peak have no influence on a It 
is also immediately evident that, if the bridging 
mechanism is such that a residual stress simply 
translates the increasing portion of the curve p(u) 
along the stress axis by o° without changing its 
shape, the matrix cracking stress must increase by 

It is useful to contrast this representation for 
the steady-state matrix cracking stress with an 
analogous interpretation for the steady state 
toughness that results when a c exceeds the peak 
value of p(u) and fibers break in the wake of the 
crack. In this case an expression for the steady 
state toughness increment bj< due to the bridging 
zone that remains over a limited region behind the 



crack tip has also been derived using the /-in- 
tegral (Budiansky et al.. 1986; Rose. 1987): 



A/ c = 2/""V(«)du 



(29) 

where u d is the crack opening displacement above 
which the ligaments no longer restrain the crack. 
Inerefore. the increase in fracture toughness is 
represented by the area beneath the curve p(u) as 
dep.cted in Fig. 4(b). Any change that increases 
this area, including modification of the peak or 
the unloading portions of the curve p(u). leads to 
an mcrease in fracture toughness. In contrast to 
the result for the steady-state matrix cracking 
stress, a reduction in stiffness of the loading por- 
tion of the curve p( U ) would usually lead to an 
increased toughness increment (provided the peak 
value of p(u), i.e. the ligament strength, is not 
decreased). The influence of residual stress on the 
toughening increment is very dependent on the 
nature of the bridging ligaments because residual 
stresses would generally be expected to affect the 
peak of the curve p(u) as well as translating the 

S ^ thC StresS ( Maf shall and Evans, 
1*88). Thus, residual stresses of a given sign can 
either increase or decrease the degree of toughen- 
ing depending on the bridging mechanisms. 
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